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A new paradigm for solar coronal heating 
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Abstract. - The solar coronal heating problem refers to the question why the temperature of the 
Sun's corona is more than two orders of magnitude higher than that of its surface. Almost 70 years 
after the discovery, this puzzle is still one of the major challenges in astrophysics. The current 
basic paradigm of coronal heating is unable to explain all the observational features of the heating. 
Here we argue that a new paradigm is required to solve the puzzle in a self-consistent manner. 
The alternative approach is based on the kinetic theory of drift waves. We show, with qualitative 
and quantitative arguments, that the drift waves have the potential to satisfy all coronal heating 
requirements. 



A self-consistent coronal heating model must fulfil a lot 
of requirements imposed by observational facts. First of 
all, it should be consistent with the measured energy losses 
in the solar corona due to conduction and radiation, i.e. 
it should a) not only provide the right amount of en- 
ergy but b) do so at the right times scales, e.g. ~ 10~ 4 
J/(m 3 s) in active regions. Moreover, it should c) include 
the source of the required energy [1], and d) work every- 
where in the corona, i.e. for all different magnetic struc- 
tures (with different heating requirements). Furthermore, 
it should be able e) to explain the observed temperature 
anisotropy [2,3] (Tj_ > Xj|), f) be more effective on ions 
than on electrons (Tj > T e ), and g) heat heavier ions more 
efficiently than lighter ions [4] . None of the proposed heat- 
ing mechanisms so far even claimed to fulfil all these model 
requirements. 

The current paradigm of coronal heating states that the 
required energy source is provided by the plasma flows 
below the solar surface and that this energy is transferred 
to the corona through the motions of the 'foot points' of 
the coronal magnetic field lines which are 'anchored' in 
this zone. Depending on the ratio r of the time scale of 
these 'driving' motions to the dynamic (Alfven) time scale, 
the current models are classified as 'wave heating' (r < 1) 
or 'magnetic reconnection (or nanoflare)' (r > 1) models 
[5]. The main challenge, however, is to explain how the 
energy is dissipated in the highly conductive corona with a 
Lundquist number (the ratio of the dissipation time to the 



Alfven time scale) around 10 13 . All proposed mechanisms 
either have a problem with the energy transport to the 
corona (the observed wave fluxes are too small) or with the 
dissipation (too little and/or too slow, or only sporadic) 
[5]. 

Most current models rely on the continuum or fluid ap- 
proximation (Magnetohydrodynamics, MHD). However, 
these models cannot really explain coronal heating com- 
pletely because i) it is clear that the actual heating takes 
place at length scales much smaller than those on which 
the (macroscopic) MHD model is justified; and ii) it is 
obvious that the observed discrepancy between ion and 
electron temperatures in the corona, as well as iii) the ob- 
served large temperature anisotropy in the inner corona 
(Tj_ > Tii), and iv) the observed preferential heating of 
the heavier ions [2] are beyond the (single!) fluid model. 

Here, we make a starting step towards the formulation 
of a new paradigm, based on the kinetic theory of the drift 
waves driven by density gradients that are omnipresent in 
the solar corona. It implies that the direct energy supply 
for the heating comes from the corona itself (from the den- 
sity gradients), though still maintained and replenished 
by some mechanisms below the surface. Those include a 
continuous restructuring of the magnetic field, implying 
consequent similar changes of the plasma density (due to 
the frozen- in conditions), and also the observed inflow of 
the plasma along the magnetic loops [6] . To some extent 
this looks similar to the currently accepted scenarios men- 
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tioned above, where the magnetic field plays an essential 
role and is assumed as given. However, the suggested new 
approach not only enables to describe these drift waves 
(which are missing in the MHD picture), but also their 
dissipation is easy to explain in the self-consistent kinetic 
model that works on the (very small) length scales at 
which the actual dissipation takes place. Actually, two 
mechanisms of energy exchange and heating will be shown 
to take place simultaneously, one due to the Landau effect 
in the direction parallel to the magnetic field, and an- 
other one, stochastic heating, in the perpendicular direc- 
tion. Moreover, this stochastic drift wave heating mech- 
anism seems to satisfy all seven above-mentioned model 
requirements. This will be proved below using only estab- 
lished basic theory, verified experimentally in laboratory 
plasmas. 

The drift mode is the only mode that is not only able to 
survive the drastically different extremes in various parts 
of the solar atmosphere but, in fact, even manages to ben- 
efit (i.e. to grow) from each of them. The driving mecha- 
nism, however, is always the same, viz. the density gradi- 
ents perpendicular to the ambient magnetic field vector. 
Numerous observations confirm the omnipresence of such 
irregularities of the plasma density across magnetic flux 
surfaces (see Fig. f ). Extremely fine density filaments and 
threads have been observed in the solar atmosphere for a 
long time now, even from ground-based observations [7] 
like those during the eclipse in f 99f , showing a slow radial 
enlargement of the structures. Contour maps [8] reveal the 
existence of numerous structures of various sizes. Filamen- 
tary structures with length scales of the order of 1 km have 
been discussed [9]. Recent Hinode observations [fO] con- 
firm that the solar atmosphere is a highly structured and 
very inhomogeneous system, with radially spreading den- 
sity filaments of various size pervading the whole domain. 
A very recent three-dimensional analysis of coronal loops 
reveals short-scale density irregularities within each loop 
separately [11]. Such density irregularities 3.1" G , 3.S 3 rule, 
associated with the magnetic field, thus creating a perfect 
environment for drift waves. The characteristic dimen- 
sions of the observed density irregularities are limited by 
the resolution of the current instruments that is presently 
a fraction of an arcsec. However, even extremely short, 
meter-size length scales can not be excluded, especially 
in the corona [12]. Hence, in dealing with drift waves, 
we may operate with density inhomogeneity length scales 
that have any value from one meter up to thousands of 
kilometers (in the case of coronal plumes). Nevertheless, 
the role of the drift wave in the problem of solar coronal 
heating has been overlooked so far in the literature, prob- 
ably due to the fact that it simply does not exist in the 
widely used MHD model. 

Regarding the wave heating analysis, in the case of the 
practically collision-less corona, an efficient mechanism 
for the transfer of energy from the wave to the plasma 
is needed. Within the framework of the drift wave ki- 
netic theory, this process develops as follows: the inter- 




Fig. 1: An image of the Sun taken by the Transition Region 
and Coronal Explorer (TRACE). 



action between the wave and the electrons is destabiliz- 
ing and the mode grows due to a Cherenkov-type inter- 
action. At the same time, however, its energy is ab- 
sorbed by the ions due to Landau damping. This may 
be seen from Ref. [13], where the drift wave properties 
within the limits k z v T i C w C k z v Te , to <C Qi, and 
Iky/kzKTe/Ti) 1 / 2 pi/ 'L n > 1, are described by the fre- 
quency 
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^di = 

— (v^, e /fl e )e z x S7±n /n , and I is the modified Bessel 
function of the first kind and of the order 0. 

Eq. (1) reveals the energy source already in the real part 
of the frequency uj r oc V±uq, while details of its growth 
due to the same source are described by Eq. (2). It is seen 
that, for a wave frequency below w* e , the growth rate can 
be written as u>i = \j e i \ — \jion\- The energy exchange 
between the plasma particles and the wave is provided by 
two mechanisms simultaneously, viz. one in the parallel 
and one in perpendicular direction. The term \~fion\ is re- 
sponsible for the Landau dissipation of the wave energy 
and, consequently, for the parallel heating of the plasma. 
So, as long as the density gradient is present, there is a 
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Fig. 2: The growth rate (2) normalized to u> r in terms of the 
parallel and perpendicular wave-lengths. The parameters are 
B = 1(T 2 T, n = 10 15 m" 3 , L n = [{dn / 'dx) / no]" 1 = 100 m, 
Ti = 10 6 K. 

continuous precipitation of energy from the wave to the 
plasma. On the other hand, the term \^ e i\ results in the 
growth of the wave, and this implies another (stochas- 
tic) heating mechanism that also involves single particle 
interactions with the wave. Remark that this process is 
described and even experimentally verified [14, 15] for a 
hot, fully ionized plasma, viz. in a tokamak. For drift 
wave perturbations of the form <f)(x) cos(k y y + k z z — u)t), 
with \k y \ 3> |fc z |, one finds the ion particle trajectory in 
the wave field from the following set of equations: 

dx/dr = T, x = k y x, T = k y y, r = VL t t, (3) 

d 2 T/dr 2 = -T + [mikl4>/ (eB 2 )] sin (T - rw/Sk) ■ (4) 

It has been shown [14] that stochastic heating takes place 
for a large enough wave amplitude, more precisely for 
a = k 2 p 2 e<j) I (kTj) > 1. The maximum achieved bulk ion 
velocity, proportional to the wave amplitude, is given by 

Vmax ^ [ky p 2 ecj)/ (nTi) + 1.9]fi i( /fc y . (5) 

Ideally, this all requires \-f e i\ > \ji on \ so that the wave 
amplitude may grow and at some moment both heating 
mechanisms may take place simultaneously. In the solar 
corona this condition can easily and quickly be satisfied 
because of the almost unlimited range of the parallel wave 
number k z , so that the ion Landau damping can be made 
small, i.e. \uj/k z \ 3> v T i, and because of the large growth 
rate as shown in the example in Fig. 2. 

In the stochastic heating process due to the drift wave, 
the ions move in the perpendicular direction to large dis- 
tances and feel the time-varying field of the wave due to 
the polarization drift v p = (dE / 'dt) / '(fijiJo) , and as a re- 
sult their motion becomes stochastic. The polarization 




Fig. 3: The temperature obtained by stochastic heating, in 
terms of the ion mass u = rrn/m v and the perpendicular wave- 
length. 

drift is in the direction of the wave vector, which em- 
phasizes the crucial electrostatic nature of the wave in 
the given heating process. Also, this stochastic heating 
is highly anisotropic, and it takes place mainly in the di- 
rection normal to the magnetic field Bq (both the x- and 
y-direction velocities are stochastic). The perpendicular 
heating in the experiment [14, 15] was larger by about a 
factor 3 compared to the parallel one, and it is exception- 
ally fast (see below). At the same time, in view of the 
mass difference and the physical picture given above, this 
heating scenario predominantly acts on the ions. 

In application to coronal magnetic structures, the indi- 
cation or 'proof that the heating really takes place would 
be: i) an ion temperature anisotropy Tjj_ S> Ti Z , ii) a pos- 
sibly higher ion temperature in comparison to electrons, 
and iii) a better heating of heavier ions. Observations 
show that i) may be taken rather as a rule than as an ex- 
ception [2,3], i.e., the perpendicular stochastic heating is 
more dominant compared to the parallel heating. There 
are also numerous indications that confirm the features ii) 
and iii). As an example we refer to graphs from Ref. [4], 
where T e < Th < Tjj e throughout the corona and the 
solar wind. 

An easy way to demonstrate that the heating can take 
place at various inhomogeneity scale-lengths L n (in other 
words in various magnetic field or density structures in 
the solar corona), is to keep the ratio X z /L n fixed. For 
example, setting L n = s X 100 m, and X z — s x 40000 m, 
where s takes values, e.g., between 1 and 1000, we calcu- 
late the frequency and the growth rate for the fixed value 
X y — 0.5 m, and we find out that the ratio Lo r /uji ~ 1. 
Note that such a variation of s may also be used to de- 
scribe the natural change of the radial density gradient 
with the increased altitude, in other words the heating 
occurs everywhere along a given flux tube. 
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Assuming an initial perturbation of the order 
e(j>/(nTi) ~ 0.01, (i.e., <fi = 0.86 V) for the parameters 
from Fig. 2 (and for ~ uj r ~ 2.5 • 10 2 Hz, A y ~ 0.5 m) 
the growth time r g = In 100/a^ till it becomes of the order 
of unity (cf> = 86 V) is about 0.02 s. For the temperature 
increased by 10 6 K (see Table 1) this implies a heating rate 
of ions of the order of 5 • 10 7 K/s, which is similar to the 
heating rate obtained in the experiments [14, 15]. Observe 
also that the magnitude of the electric field which we are 
dealing with is of the same order as in those experiments. 

The stronger heating of heavier ions (see Fig. 3) can 
be understood from Eq. (5) and after expressing the ef- 
fective temperature in terms of the ion mass T e ff(rrii) = 
m i v max/(3 K )- From the derivative dT e j / (mj) / dnii > 0, 
it follows that the heating increases with the ion mass 
if k*pf(e<f)/KTi) 2 > 1.9. For = 60 V we have the 
normalized temperature Ti(A y , p) = 0.881 + 0.057/x/A 2 , + 
1.78A 2 / M . 

For the same parameters as above, the maximum en- 
ergy released per unit volume is T, max = nomiV^ nax / 2 = 
0.04 J/m 3 . The energy release rate T max = Y< max /T g ~ 
1 J/(m 3 s) amounts to 4 orders of magnitude above the 
necessary value. However, for L n — 100 km (i.e., set- 
ting s = 1000) we obtain lu, = 0.13 Hz, uj r = 0.254 Hz, 
T g = 32.6 s, and consequently T max = 1.2 • 10~ 3 J/(m 3 s), 
that is about one order of magnitude above the heating 
rate accepted as necessary. Similar estimates may be done 
for still larger L n , yet the conditions under which the pre- 
vious expressions are derived become violated and a nu- 
merical approach is required in this case. In reality the 
collisions and nonlincarity lead to the flattening of the 
density profile in the region occupied by the wave [16], 
that should result in the saturation of the growth. Also, 
energy diffusion in the perpendicular direction should act 
in the same way. In addition, the drift wave is as a rule 
coupled to the Alfven wave [17], with the coupling propor- 
tional to kyPi. All these effects will more effectively act on 
short scales, and the actual values for T are expected to 
be below T max . Therefore, the apparently too big release 
of energy at short scales, as formally obtained above, may 
in reality be considerably reduced. Clearly, more accurate 
estimates and more detailed description may be obtained 
only numerically. 

To summarize, the proposed mechanism is based on 
a novel paradigm that allows a self-consistent solution 
model. The heating mechanism implies instabilities on 
time and spatial scales that are currently not directly ob- 
servable by space probes. However, all the effects pre- 
sented here are directly experimentally verified under lab- 
oratory conditions. Their indirect confirmation in the con- 
text of the solar corona seems to be also beyond doubts. 
This is because the consequences of the heating process, as 
enlisted earlier in the text (temperature anisotropy, bet- 
ter heating of heavier ions, hotter ions than electrons), are 
indeed verified by satellite observations. 



Table 1: Plasma heating for hydrogen ions for two perpendic- 
ular wavelengths and for two values of the wave amplitude 
(j>. The values in brackets are for helium. The maximum 
stochastic velocity is given by Eq. (5) and the effective tem- 
perature obtained by heating is mi„ / '(3k) . We used the 
same starting set of parameters as earlier, viz. L n — 100 m, 
\ z = 20 km, n = 10 15 m~ 3 , T = 10 6 K. For shorter wave- 
lengths the heating is stronger for helium because the condition 
ky-piie^/nTi) 2 > 1.9 (see Fig. 3) is easily satisfied. 



<t> = 60 [V] 


<t> = 80 [V] 




A, [m] T eff [K] 


Teff [K] 




0.5 1.56 x 10 b (1.91 x 10°) 


2.03 x 10° (2.92 x 


10°) 


1 2.72 x 10 6 (1.56 x 10 6 ) 


3.06 x 10 6 (2.03 x 


10 6 ) 



* * * 

These results are obtained in the framework of 
the projects G.0304.07 (FWO-Vlaanderen) , C 90205 
(Prodex 9), and GOA/2009-009 (K.U.Leuven). TRACE 
(Fig. 1) is a mission of the Stanford-Lockheed Institute 
for Space Research, and part of the NASA Small Explorer 
program. 
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